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1.  Introduction 


Designing  components  fabricated  from  ceramic  materials  requires  the  capability  to  extrapolate 
failure  statistics  from  test  specimens  to  the  components,  i.e.,  one  must  be  able  to  account  for  size 
effects  commonly  observed  in  ceramic  material  systems.  The  component  of  interest  in  this 
endeavor  is  a  ceramic  gun  barrel.  Before  the  program  began,  planning  that  supported  this  effort 
anticipated  the  need  to  capture  defects  spatially  distributed  around  the  circumference  of  the 
barrel,  as  well  as  defects  spatially  distributed  down  the  length  of  the  barrel.  This  focused 
attention  on  two  types  of  test  specimen  geometries.  Specifically,  a  C-ring  geometry  was  adopted 
to  investigate  circumferential  flaws,  and  a  sectored  flexure  specimen  subjected  to  four-point 
bending  will  be  used  to  investigate  longitudinally  distributed  flaws.  A  review  of  currently 
available  analytical  methods  is  presented  for  the  C-ring  specimen,  and  results  from  extensive 
computer  modeling  efforts  are  presented  herein  for  both  test  specimen  geometries.  The  test 
specimen  geometries  were  defined  by  the  U.S.  Army  Research  Laboratory  (ARL),  and  the  intent 
was  to  detennine  the  effective  volume  and  effective  area  (respectively  referred  to  as  ky‘V and 
Jca'A  in  the  ceramics  literature)  as  a  function  of  the  Weibull  modulus  and  geometry.  Here,  the  ky 
and  kA  are  stress  gradient  factors,  and  V  and  A  are  gauge  volume  and  gauge  area,  respectively. 
Calculation  of  kpV  and  k,\ -A  required  numerical  computations  using  both  finite  element  analysis 
and  Weibull  analysis.  The  ANSYS1  finite  element  program  was  used  to  generate  stress  results 
for  this  effort,  while  the  National  Aeronautics  and  Space  Administration  Glenn  Research  Center 
Ceramics  Analysis  and  Reliability  Evaluation  of  Structures  (CARES)/Life  algorithm  was  used  to 
conduct  Weibull  analyses. 


2.  Effective  Volumes  and  Areas  Using  Finite  Element  Analysis 


If  the  random  variable  representing  uniaxial  tensile  strength  of  an  advanced  ceramic  is 
characterized  by  a  two-parameter  Weibull  distribution,  then  the  failure  probability  for  a  test 
specimen  fabricated  from  this  advanced  ceramic  is  given  by  the  cumulative  distribution  function 


Pf  =  1  -exp 


V  <Je  ) 


>0 


(1) 


Pf  =  0  <7  <0 

1J  max  — 


(2) 


'ANSYS®,  which  is  not  an  acronym,  is  a  registered  trademark  of  ANSYS,  Inc. 
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in  which 

Pf  =  the  probability  of  failure 

(Tmax  =  maximum  uniaxial  tensile  stress  in  a  test  specimen  at  failure 

<Jq  =  the  Weibull  characteristic  strength 

m  =  Weibull  modulus 


Note  that  the  Weibull  characteristic  strength  depends  on  the  geometry  of  the  test  specimen.  In 
addition,  the  Weibull  characteristic  strength  has  units  of  stress.  Strength-controlling  flaws  can  be 
spatially  distributed  over  the  volume,  the  surface  area,  or  edge  length  (three-,  two-,  or  one¬ 
dimensional  entities,  respectively)  of  a  test  specimen.  If  the  strength-controlling  flaws  are 
spatially  distributed  over  the  volume  of  the  test  specimen,  the  characteristic  strength  typically  is 
designated  as  (<Je)v,  and  the  Weibull  modulus  is  typically  designated  as  mv.  If  the  strength¬ 
controlling  flaws  are  spatially  distributed  over  the  area  of  the  test  specimen,  the  characteristic 
strength  is  designated  as  (<Jq)a,  and  the  Weibull  modulus  is  designated  as  mA.  If  the  strength¬ 
controlling  flaws  are  distributed  over  the  edge  of  a  test  specimen,  then  the  characteristic  strength 
is  designated  as  (<Jo)l,  and  the  Weibull  modulus  is  designated  as  /«/_.  It  has  been  recognized  for 
several  decades  that  edge-based  strength-limiting  flaws  in  ceramics  can  be  particularly 
debilitating  to  failure  stress.  Typically,  edges  in  ceramic  specimens  and  components  are 
chamfered  or  rounded  for  that  reason  (i.e.,  surface-  and  volume-based  flaws  are  controlling); 
however,  edge  alteration  is  not  always  practical  in  some  applications,  so  the  probability  of  failure 
attributable  to  these  defects  must  be  accounted  for.  Obviously,  for  tubular  or  rounded  specimens 
and  components,  the  consideration  of  edge-based  strength-limiting  flaws  is  not  necessary. 


An  alternate  expression  for  the  probability  of  failure  is  given  by 


Pf  =  1 -ex  p 


1 


dV 


<J>  0 


(3) 


=  Pf  =  0  <J  <0 


(4) 


The  integration  within  the  exponential  function  is  performed  over  all  tensile  regions  of  the 
specimen  volume  if  the  strength-controlling  flaws  are  randomly  distributed  through  the  volume 
of  the  material.  The  integration  is  sometimes  performed  over  an  effective  gauge  section  instead 
of  over  the  total  volume  of  the  test  specimen.  In  equation  3,  my  is  the  Weibull  modulus 
associated  with  strength-controlling  flaws  distributed  through  the  volumes  of  the  specimens. 

Also  in  equation  3,  (<Jq)v  is  the  Weibull  material  scale  parameter  and  can  be  described  as  the 
Weibull  characteristic  strength  of  a  specimen  with  unit  volume  loaded  in  unifonn  uniaxial 
tension.  The  Weibull  material  scale  parameter  has  units  of  stress(volume)7/"v  when  the  strength¬ 
controlling  flaws  are  distributed  through  the  volume  of  the  material. 
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At  this  point,  we  assume  that  the  Weibull  modulus  can  be  obtained  from  the  failure  data.  Thus, 
in  equation  3  we  take 

mv  =  mv  (5) 


in  which  the  tilde  signifies  that  the  value  of  my  is  an  estimate  obtained  from  failure  data,  and  we 
identify 


as  a  parameter;  then 


Rearranging  yields 


Now  let 


(<Jo)v  ~  <?* 


(6) 

(V) 

(8) 

(9) 


in  which  V  can  be  identified  as  the  specimen  gauge  volume.  Then 


pf  =  1  -  exp 


(a  ^ 

mv 

-kv-V 

max 

V  cr*  y 

(10) 


Finally,  solve  equation  10  for  the  effective  volume,  ky'V.  For  the  time  being,  treat  a*  as  a 
parameter. 

■Mi  —  pf) 


K  -V  =  - 


\mv 


(11) 


V  o*  ) 


At  this  point,  we  have  outlined  the  theory  that  supports  an  approach  that  could  use  finite  element 
analysis  and  Weibull  analysis  to  evaluate  the  effective  volume,  ky'V,  for  a  given  specific  test 
specimen  geometry.  The  following  steps  are  necessary  to  evaluate  ky'V: 

1 .  Estimate  mv  from  failure  data.  In  the  case  of  the  parametric  study  presented  later,  a  range 
of  values  was  assumed  for  mv . 


2.  Use  finite  element  analysis  and  Weibull  analysis  to  numerically  evaluate  the  following 
expression  for  a  given  specimen  geometry  from  equation  7: 


Pf 


1  -exp 


-1 


(7) 
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Here,  mv  from  step  1  is  used  with  an  arbitrarily  selected  value  of  o*.  Weibull  analysis  will 
numerically  determine  Pt  for  the  (arbitrary)  magnitude  of  the  load  applied  to  the  specimen. 


3.  With  a*,  Pf,  and  <Jmax  (obtained  from  the  finite  element  analysis),  solve  the  following 
expression  from  equation  1 1 : 


ky-V  = 


-h(l  -  Pf) 


\my 

) 


(11) 


Note  that  an  increase  in  the  applied  load  will  increase  both  Pf  and  <7max  so  that  ky'V  remains 
constant  for  a  given  specimen  geometry.  In  addition,  changing  a*  will  affect  the  value  of  Pf 
obtained  from  step  2  so  that  ky'V  remains  constant  for  a  given  material  characteristic  strength. 
The  only  parameter  other  than  the  geometry  that  will  affect  the  value  of  ky'V  is  the  Weibull 
modulus  rhv . 


3.  Finite  Element  Analysis  of  the  C-Ring  Specimen 


In  this  section  and  several  sections  that  follow,  the  methods  of  calculating  the  effective  volume 
and  area  from  finite  element  modeling  for  the  C-ring  specimen  geometry  are  compared  with 
calculations  obtained  with  the  published  closed  form  solutions  provided  by  Jadaan  et  al.  (1991) 
and  presented  in  American  Society  for  Testing  Materials  (ASTM)  C-1323.  Throughout  this 
report,  the  solution  in  Jadaan  et  al.  (1991)  will  be  referred  to  as  the  closed  form  solution.  For  this 
project,  the  C-ring  geometry  ideally  has  an  outside  diameter  (OD)  of  33  mm  and  an  inside 
diameter  (ID)  of  24  mm.  The  specimen  length  for  the  C-ring  was  8  mm,  but  several  other 
lengths  (specifically  2  mm  and  4  mm)  were  examined  to  indicate  certain  analytical  trends.  Finite 
element  analysis  provides  a  numerical  approximation  of  the  stress  field  and  the  geometry  of  a 
component  being  analyzed.  Since  the  probability  of  failure  depends  on  stress  and  geometry,  a 
mesh  sensitivity  study  was  conducted  for  the  C-ring  geometry  provided.  Five  meshes  were 
analyzed.  These  meshes  and  the  corresponding  number  of  elements  are  identified  as  follows: 

1.  Most  coarse  mesh:  864  elements 

2.  Coarse  mesh:  2,800  elements 

3.  Coarse  to  medium  mesh:  6,860  elements 

4.  Medium  mesh:  16,128  elements 

5.  Fine  mesh:  19,872  elements 

The  most  coarse  mesh  and  the  most  refined  mesh  appear  in  figures  1  and  2. 
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ANSYS  SOLID952  elements  were  used  to  model  the  C-ring.  A  series  of  simulations  was 
conducted  on  the  ARL  specimen  geometry  with  the  meshes  cited  previously.  Force  and 
displacement  boundary  conditions  were  imposed  at  the  top  of  the  C-ring.  For  a  given  specimen 
geometry,  both  approaches  yielded  identical  results.  In  this  report,  displacement  models  are 
presented. 


Figure  1.  Most  coarse  C-ring  mesh. 


Figure  2.  Most  refined  C-ring  mesh. 


7  ....  . 

“SOLID95  is  a  higher  order  version  of  the  3-D  eight-node  solid  element  SOL1D45.  It  can  tolerate  irregular 
shapes  without  as  much  loss  of  accuracy. 
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The  first  principal  stress  contour  plot  for  the  ARL  specimen  geometry  is  depicted  in  figure  3. 

This  contour  plot  is  representative  of  the  five  meshes  cited  before.  Although  the  magnitude  of 
the  first  principal  stress  is  not  of  consequence  to  this  discussion,  the  distribution  of  the  stresses  is. 
Note  the  distinct  parabolic  distribution  of  stress  along  the  length  of  the  C-ring  in  figure  3. 


Figure  3.  Principal  stress  contour  map. 


From  the  closed  form  solution,  the  maximum  tensile  stress  (which  occurs  at  the  outer  surface  of 
the  bend  face)  of  a  C-ring  with  an  applied  compressive  load  is  given  in  this  article  as 


in  which 

P  =  load 


P  R 

<7  = - 

max  *  . 

bt  r 


o  \  a 


R 


(12) 


t  =  radial  thickness 

b  =  length 

rG  =  outside  radius 
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ra  =  average  radius 
r,  =  inner  radius 


and 


R 


In 


V  ri  J 


A  schematic  drawing  depicting  the  various  geometries  identified  is  provided  in  figure  4. 


Figure  4.  Schematic  of  C-ring  specimen  geometry. 


(13) 


An  elasticity  solution  is  also  available  (see  Timoshenko  and  Goodier,  1970).  The  stress  equation 
obtained  from  the  elasticity  solution  has  never  been  used  to  compute  effective  volumes  or  areas 
for  the  C-ring  geometry.  The  elasticity  equations  are  presented  here  for  verification  purposes 
only. 

The  strength  of  materials  solution  and  the  elasticity  solution  are  two-dimensional  solutions,  i.e., 
the  solutions  are  correct  in  the  limit  as  the  length  of  the  C-ring  approaches  zero.  Jadaan  et  al. 
(1991)  assumed  that  this  solution  was  valid  through  the  length  of  the  C-ring  in  order  to  derive 
closed  fonn  solutions  for  ky-  V  and  kA  -A .  This  assumption  is  valid  during  limiting  geometric 
constraints.  Figure  5  depicts  the  variation  of  the  maximum  stress  along  the  length  of  several  C- 
rings  for  a  number  of  values  of  b/(rG  -  r,). 
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Figure  5.  Maximum  circumferential  tensile  stress  as  a  function  of  specimen  geometry. 


The  trend  to  a  uniform  stress  distribution  is  readily  apparent  in  figure  6.  Here,  the  stress  results 
are  presented  as  a  function  of  normalized  (with  respect  to  b)  position  along  the  length  of  the  bend 
face. 


Figure  6.  Maximum  circumferential  tensile  stress  as  a  function  of  specimen  geometry  normalized  with 
respect  to  specimen  length. 
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Clearly,  the  uniform  stress  distribution  is  not  satisfied  until 

f  r  A 

<  0.5 


\ro~ri  J 


(14) 


This  result  has  practical  implications  for  the  testing  community.  The  ARL  C-ring  geometry  was 
designed  to  comply  with  geometric  constraints  specified  in  ASTM  C-1323  in  which 


1  < 


Kro~ri  J 


(15) 


The  ARL  C-ring  geometry  has  a  length-to-thickness  ratio  of  1.78  (in  which  b  =  8  mm  was 
specified  and  the  inner  and  outer  diameters  could  not  be  adjusted),  which  falls  within  the 
recommended  limit  suggested  by  the  ASTM  standard.  From  the  stress  results  presented  in 
figure  4,  we  feel  that  the  ASTM  C-28  committee  should  revisit  this  recommendation.  If  the 
inequality  cited  in  equation  1 5  has  practical  testing  implications  (specimen  alignment,  stability, 
etc.),  then  a  caveat  should  appear  within  the  standard  that  states  the  expressions  for  ky’V and  kA-A 
in  the  standard  are  in  appreciable  error. 


4.  Finite  Element  Analysis  of  the  Sectored  Flexure  Specimen 


The  sectored  flexure  bar  geometry  was  based  on  a  tube  with  dimensions  of  33  mm  OD,  24  mm 
ID,  and  100  nun  length.  The  four-point  bend  test  geometry  consisted  of  an  outer  support  span  of 
80  mm  and  an  inner  load  span  of  40  mm.  The  cross  section  of  the  contoured  specimen  was  a  45- 
degree  sector  machined  from  the  tube  geometry  (see  figure  7). 


Figure  7.  Sectored  tube  geometry. 
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A  flat  surface  was  machined  on  the  compressive  side  of  the  contoured  flexure  specimen, 
producing  a  specimen  with  a  maximum  thickness  of  4  mm.  Bottom  support  fixtures  (see 
figure  8)  were  assumed  to  have  of  an  “hourglass”  shape  with  the  same  radius  of  curvature  as  the 
contoured  flexure  specimen,  promoting  line  loading  which  was  modeled  via  finite  element 
analysis. 


Figure  8.  One  of  two  contoured  support  fixtures  (courtesy  of 
S.  Goodrich  at  the  University  of  Dayton 
Research  Institute). 

Initially,  models  that  have  force  and  displacement  boundary  conditions  imposed  at  the  inner  load 
supports  were  considered.  The  inner  load  supports  are  situated  along  the  top  of  the  specimen  40 
mm  apart  (as  depicted  in  figure  7).  The  force  boundary  conditions  (outlined  in  red)  are  depicted 
in  figure  9,  which  represents  a  model  of  the  entire  test  specimen  geometry.  The  hourglass 
supports  for  the  bottom  of  the  specimen  are  depicted  in  light  blue. 

For  the  applied  displacement  model,  a  displacement  of  0. 1  mm  was  imposed  on  all  nodes  across 
the  top  of  the  model  at  the  inner  supports.  Unfortunately,  unwarranted  stresses  were  induced 
when  this  unifonn  displacement  was  applied. 

The  stresses  corresponding  to  one  of  these  lines  of  displaced  nodes  are  depicted  in  figure  10. 
Nodes  near  the  centerline  (relative  to  x-axis)  of  the  specimen  exhibited  compressive  av  stresses 
(denoted  as  SY  in  the  figure  and  by  ANSYS)  which  were  expected  because  of  the  compressive 
nature  of  the  loading  at  these  nodes.  However,  the  <Jy  stresses  became  tensile  near  the  edges  of 
the  cross  section.  This  indicates  that  these  near-edge  nodes  sought  to  displace  farther  than  the 
applied  displacement  of  0.1  mm  and  were  actually  “held  back”  by  the  proscribed  displacement, 
thereby  causing  tensile  stresses.  Basically,  the  model’s  top  surface  has  a  tendency  to  curl  down 
at  the  edges.  The  geometry  associated  with  this  curl  cannot  be  calculated  a  priori.  This  strongly 
points  to  the  use  of  contact  stress  elements  in  order  to  model  the  test  specimen  during 
displacement  conditions.  The  use  of  contact  stress  elements  was  outside  the  present  work  scope, 
so  attention  was  focused  on  modeling  the  test  specimen  by  force  boundary  conditions. 
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Figure  9.  Full  model  with  applied  forces. 


To  ensure  that  the  force  boundary  conditions  model  was  behaving  as  expected,  displacements 
were  monitored  at  several  locations  in  the  model.  These  displacements  were  found  to  be  con¬ 
sistent  with  a  component  undergoing  four-point  bending,  and  the  decision  was  made  to  ascertain 
ky'V  and  kA  -A  via  applied  force  models.  Thus,  the  results  presented  in  the  remainder  of  the  report 
represent  finite  element  models  with  applied  force  boundary  conditions. 


Once  the  boundary  conditions  for  the  force  model  were  established  and  verified,  the  stress  state 
obtained  from  the  finite  element  analysis  was  verified.  A  strength-of-materials  solution  was 
adopted  for  the  closed  form  solution  of  the  maximum  tensile  stress  (which  occurs  at  the  outer 
surface  of  the  bend  face)  in  order  to  make  this  comparison.  The  closed  form  stress  is  given  by 
the  expression 


o 


max 


Me 

I 


(16) 
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Nodal  Stress  Profile 
Load  Line  (Y=4,  Z=30) 
Applied  Displacement  UY  =  -0.1 


Width  (mm) 

Cross  Section  Centerline  @  X  =  0 

Figure  10.  Nodal  stresses  along  load  support  (Z  =  20  mm  from  center)  across  top  surface  (Y  =  4  mm  from 
bottom):  displacement  loading. 

Additionally,  an  independent  calculation  for  the  moment  of  inertia  was  compared  to  the  value 
obtained  from  ANSYS.  For  a  1000-N  load  (500  N  applied  to  each  of  the  load  lines),  the  closed 
form  solution  yields  a  maximum  tensile  stress  of  445.66  MPa.  The  maximum  stress  obtained 
from  the  most  refined  mesh  (see  discussion)  with  ANSYS  was  445.76  MPa.  Note  that  the 
material  response  was  assumed  to  be  elastic  throughout.  Moreover,  the  authors  wish  to  point  out 
that  it  is  not  germane  to  the  analysis  of  ky-  V  and  kA'A  that  the  loads  and  corresponding  stresses 
are  too  high  relative  to  the  candidate  materials  undergoing  consideration  in  the  gun  barrel 
program.  For  a  given  Weibull  modulus,  the  parameters  kyV  and  kA-A  depend  only  on  load 
configuration  and  test  specimen  geometry,  not  on  ultimate  strength  or  elastic  material  properties. 

In  preliminary  analyses,  results  for  ky-  V  and  kA -A  were  generated  for  the  sectored  flexure  bar 
modeled  in  its  entirety,  i.e.,  the  full  model.  In  later  analyses,  the  models  were  simplified  via  half 
and  quarter  symmetry  in  order  to  save  on  computational  resources.  For  identical  mesh  densities, 
the  full,  half,  and  quarter  models  produced  identical  kyV  and  kA-A  parameters.  Since  the  quarter 
symmetry  models  produced  the  least  demand  on  computational  resources,  this  model  was  used 
for  all  subsequent  analyses. 
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Finite  element  analysis  provides  a  numerical  approximation  of  the  stress  field  and  the  geometry 
of  a  component  being  analyzed.  Since  the  probability  of  failure  depends  on  stress  and  geometry 
and  stress  depends  on  mesh  density,  a  mesh  sensitivity  study  was  conducted  for  the  sectored 
flexure  bar.  Gradually  increasing  the  mesh  density  allows  for  an  asymptotic  approach  to 
estimations  of  the  ky'V  and  kt  -A  parameters. 

ANSYS  SOLID95  elements  were  used  to  model  the  sectored  flexure  bar.  A  CARES  surface 
macro  was  executed  for  each  model  to  compute  the  surface  of  the  bar  for  subsequent  use  in  the 
calculation  of  k a 'A .  Five  meshes  were  analyzed  with  quarter  symmetry.  These  meshes  and  the 
corresponding  number  of  elements  are  identified  in  table  1 . 


Table  1.  C-ring  elements  according  to  mesh  density 


Mesh  Density 

Volume  Elements 

Surface  Elements 
(CARES  generated) 

Most  Coarse 

256 

208 

Coarse 

1,500 

800 

Coarse/Medium 

4,480 

1,552 

Medium 

10,500 

2,410 

Fine 

21,000 

3,800 

The  most  coarse  mesh  and  the  fine  mesh  appear  in  figures  1 1  and  12.  The  (7Z  stress  (denoted  SZ 
by  ANSYS)  contour  plot  for  the  quarter  symmetry  mesh  is  depicted  in  figure  13.  This  stress  is 
not  a  principal  stress,  but  the  linear  distribution  of  this  stress  is  expected  in  beam  bending.  This 
contour  plot  is  representative  of  the  five  meshes  cited  before.  Although  the  magnitude  of  this 
stress  is  not  of  consequence  to  this  discussion,  the  distribution  of  the  stresses  is.  The  stress 
contour  plot  is  further  confirmation  that  the  finite  element  model  correctly  captures  the  expected 
response  of  the  test  specimen. 

Figure  14  depicts  a  contour  plot  of  the  first  principal  stress.  Since  the  CARES/Life  reliability 
analysis  will  use  tensile  principal  stresses,  this  figure  indicates  the  region  over  which  the 
reliability  integration  will  be  computed.  Note  that  the  tensile  portion  of  the  first  principal  stress 
“dies”  in  the  vicinity  of  the  outer  support  span.  Thus,  the  volume  outside  the  outer  support  will 
contribute  very  little  to  the  ky'V  or  kA -A  calculation. 
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Figure  11.  Most  coarse  quarter  symmetry  mesh  with  applied  forces. 


Figure  12.  Most  refined  quarter  symmetry  mesh  with  applied  forces. 
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Figure  13.  Z-axis  stress  contour  plot. 
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Figure  14.  Principal  stress  contour  plot. 
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5.  C-Ring  Effective  Volume 


Before  the  stress  results  cited  previously  were  established,  the  closed  form  solutions  for  kv’V and 
kA‘A  cited  in  Jadaan  et  al.  (1991)  were  considered  benchmarks  which  the  numerical  approach 
using  finite  element  analysis  and  Weibull  analysis  had  to  achieve.  The  closed  form  C-ring 
solution  (compressive  loading)  for  ky-V is 


in  which 


K-v 


(<Je) 


dV 


max  J 


(17) 


fe  = 
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b  rom  fe 

fr 

rn/2  ,  . 
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(  \ 

m 

1  'a 

1 —m  i 

r  dr 

(19) 

Vo~ra J 

As  indicated  before,  the  computation  of  kv’V  depends  on  the  Weibull  modulus.  This  dependency 
is  depicted  in  figure  15  where  the  kv’V  value  for  each  mesh  refinement  is  plotted  as  a  function  of 
the  assumed  Weibull  modulus. 

This  plot  is  cited  in  Jadaan  et  al.  (1991)  and  the  relative  proximity  of  the  curves  generated  by 
numerical  analysis  with  respect  to  the  closed  fonn  solution  is  somewhat  reassuring.  It  is 
interesting  to  note  that  the  most  coarse  mesh  tracks  the  closed  form  solution  the  best.  However, 
a  more  representative  graph  would  be  the  percentage  difference  in  the  numerical  computation  of 
kv’V  with  respect  to  the  closed  form  solution.  This  infonnation  is  presented  in  figure  16. 

In  figure  16,  it  is  quite  evident  that  the  most  coarse  mesh  is  closely  tracking  the  closed  form 
solution.  This  is  counter-intuitive  and  leads  to  examining  the  numerical  results  in  a  manner 
depicted  in  figure  17.  The  graph  in  figure  17  indicates  that  the  numerical  computations  are 
converging  on  an  asymp-totic  value  different  than  the  closed  form  solution  as  the  mesh  is 
refined.  The  results  presented  in  figure  18  correspond  to  a  Weibull  modulus  of  10,  which  is 
a  value  that  should  be  representative  of  a  material  undergoing  consideration  in  the  gun  barrel 
project. 

The  trends  cited  previously  exist  not  only  for  the  Weibull  modulus  but  also  for  the  specimen 
geometry.  Consider  the  results  depicted  in  figure  19  where  the  percentage  difference  in  the 
numerical  computation  of  kv’V  with  respect  to  the  closed  form  solution  is  graphed  as  a  function 
of  geometry. 
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Coarsest  — ■ —  Coarse  — * —  Cor/Med  — x —  Medium  — x —  Fine - Closed  Form 


Figure  15.  Effective  volume  ( kv-V)  as  a  function  of  the  Weibull  modulus  and  mesh. 


— ♦ — Coarsest  — ■ — Coarse  a  Cor/Med  x  Medium  x  Fine - Closed  Form 


Figure  16.  Percentage  difference  in  (kv-V)  as  a  function  of  mesh  size  and  Weibull  modulus. 
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Figure  17.  Effective  volume  as  a  function  of  mesh  size  for  m  =  10. 

Here,  the  percentage  difference  exhibited  for  the  most  refined  meshes  (which  represent  the  most 
accurate  numerical  computation)  increases  as  the  length  (b)  of  the  specimen  increases  (the  inner 
and  outer  radii  are  held  constant  in  this  study).  As  long  as  the  value  of  b/(ra  -  r,)  stays  below  l, 
there  is  reasonable  correlation  between  the  two  approaches. 


m  =  10 


o 

a> 

a. 

-40%  - 

0.25  0.5  0.75  1  1.25  1.5  1.75  2 

b/(rQ  -  n) 


Figure  19.  Percentage  difference  in  effective  volume  relative  to  the  closed  form  solution 
as  a  function  of  specimen  geometry  for  m  -10. 
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6.  C-Ring  Effective  Area  (kA) 


The  expression  for  kA‘A  cited  in  Jadaan  et  al.  (1991)  is  as  follows 


=  brofe  +2  r0m  fe  fr 


(20) 


in  which  fe  and  [■  are  given  by  equations  18  and  19. 


Similar  asymptotic  behavior  was  found  in  the  numerically  generated  values  of  the  effective  area 
parameter  kA -A .  This  convergence  to  a  value  different  than  the  closed  fonn  solution  is  depicted 
in  figure  19. 


FEA  — ■ —  Closed  Form 


Total  Number  of  Elements 


Figure  19.  Effective  area  (kA  -A)  as  a  function  of  mesh  size  for  m  =  10. 


A  quantitative  assessment  in  the  difference  between  the  two  values  over  a  range  of  Weibull 
modulus  values  is  given  in  figure  20. 
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m  =  10 


-♦ —  Coarsest  — ■ —  Coarse  —a —  Cor/Med  — x —  Medium  — x —  Fine - Closed  Form 


Figure  20.  Percentage  difference  in  effective  area  (kA  A)  relative  to  the  closed  form  solution  as 
a  function  of  specimen  geometry  for  m  =  10. 


7.  Sectored  Flexure  Specimen  Effective  Volume 


Closed  form  solutions  for  ky  V  and  k  \ -A  do  not  exist  for  this  specimen  geometry.  As  indicated 
before,  the  computation  of  kyV  depends  on  the  Weibull  modulus.  This  dependency  is  depicted 
in  figure  2 1  where  the  kyV  values  for  each  mesh  refinement  are  plotted  as  a  function  of  the 
assumed  Weibull  modulus.  The  graph  in  figure  22  indicates  that  the  numerical  computations  are 
converging  to  an  asymptotic  value  for  kyV,  given  a  Weibull  modulus  of  10.  This  value  of  the 
Weibull  modulus  is  representa-tive  of  the  candidate  ceramic  materials  undergoing  consideration 
in  this  program. 

The  asymptotic  trends  indicated  in  figure  22  were  not  extended  to  yield  a  single  value  for  kyV 
since  computational  resources  were  exhausted  at  21,000  elements  using  a  quarter  symmetry 
finite  element  model.  Clearly,  the  asymptotic  trend  is  there,  but  unlike  the  Ua’A  computations 
where  a  single  pronounced  value  is  evident  (see  following  discussion),  we  can  only  surmise  a 
ky-  V  value  at  this  point.  One  or  more  mesh  refinements  would  be  necessary  to  confirm  an 
asymptotic  value. 
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Coarsest  — ■ —  Coarse  —a —  Cor/Med  — x —  Medium  — * —  Fine 


Figure  21.  Effective  volume  (kv-V)  as  a  function  of  the  Weibull  modulus  and  mesh  density. 


Figure  22.  Effective  volume  as  a  function  of  mesh  density  for  m  =  10. 
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8.  Sectored  Flexure  Specimen  Effective  Area 


As  indicated  before,  the  computation  of  kA‘A  depends  on  the  Weibull  modulus.  This  dependency 
is  depicted  in  figure  23  where  the  kA'A  value  for  each  mesh  refinement  is  plotted  as  a  function  of 
the  assumed  Weibull  modulus. 


Figure  23.  Effective  area  (kA  A)  as  a  function  of  the  Weibull  modulus  and  mesh  density. 

Similar  asymptotic  behavior  was  found  in  the  numerically  generated  values  of  the  effective  area 
parameter  k a  ’A .  This  convergence  is  depicted  in  figure  24.  Here,  the  asymptotic  value  for  kf-A 
(=  196.8  mm2)  is  clearly  evident. 
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9.  Conclusions 


The  closed  form  solutions  for  ky-  V  and  kA  -A  are  not  valid  for  the  ARL  C-ring  specimen  geometry. 
With  the  current  specimen  geometry,  a  number  of  finite  element  meshes  must  be  evaluated  to 
ensure  asymptotic  convergence  of  ky'V and  kt-A  for  a  range  of  Weibull  modulus  values.  To 
obtain  convergence  with  the  closed  form  solution,  the  C-ring  must  be  made  more  narrow.  This 
may  introduce  experimental  difficulties  with  alignment  and/or  stability  during  testing. 

Consider  the  error  introduced  in  component  reliability  calculations  if  the  closed  form  solution  for 
ky’V  is  used  to  analyze  failure  data.  The  discussion  that  follows  would  also  be  pertinent  for  kA  'A . 
Given  the  trends  illustrated  in  the  previous  sections,  the  ky'V  computations  are  approaching 
theoretically  correct  values  in  an  asymptotic  fashion.  If  the  closed  fonn  ky'V  is  in  error  relative 
to  a  theoretically  correct  value,  i.e., 

(kV  -V) dosed  form  =  ^Or  ■  (ky  '  V) ^  (21) 

in  which  the  subscript  “true”  refers  to  theoretically  correct,  and  Error  is  an  error  coefficient  that 
can  be  greater  than  or  less  than  1 ,  then 

=  °AErrorr~(kv-V)]2  (22) 

This  leads  to 
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or 

(P, )  =  1  -  (l  -  (p.)  )VEmr  (24) 

x  J  'closed  form  x  x  J  'true  ' 

whereas  in  equation  22 

=  v,(K-rZ:  05) 

is  the  theoretically  correct  Weibull  material  scale  parameter.  In  tenns  of  reliability,  equation  24 
can  be  expressed  as 

K^M„=(Kj'Cr"’r  (26) 

Note  that  in  the  limit  as  the  error  coefficient  becomes  large  (i.e.,  >1),  the  closed  form  reliability 
approaches  1,  and  the  closed  form  probability  of  failure  trends  to  0.  This  trend  is  non-conserva¬ 
tive.  On  the  other  hand,  as  the  error  coefficient  approaches  0  in  the  limit  (from  1),  then  the 
reliability  and  probability  of  failure  computations  obtained  from  the  closed  form  approach 
become  overly  conservative. 

Specifically,  assume  that  the  error  coefficient  assumes  a  value  Error  =  1.25.  This  value  is 
similar  to  that  associated  with  the  ARL  C-ring  geometry  for  a  Weibull  modulus  of  10.  Next, 
assume  that  the  component  reliability  based  on  the  theoretically  correct  value  of  o0v  is  0.99999, 
or  1  failure  in  100,000.  The  expression  given  by  equation  26  would  then  yield  a  component 
reliability  of  0.999992,  or  1  failure  in  125,000.  Thus,  the  component  probability  of  failure  is  in 
error  by  nearly  20%  (0.000008  as  opposed  to  the  theoretically  correct  value  of  0.00001).  The 
reliability  is  overestimated  (and  the  probability  of  failure  is  underestimated),  which  is  an  non¬ 
conservative  error.  Thus,  using  an  incorrect  effective  volume  or  area  to  analyze  failure  data  can 
cascade  into  component  design  and  can  yield  deleterious  results. 

The  previous  discussion  must  be  tempered  by  the  fact  that  loss  in  fidelity  in  calculating 
component  reliability  can  be  caused  by  several  sources.  These  sources  include  (a)  sampling 
error  associated  with  estimating  Weibull  parameter  from  failure  data,  (b)  experimental  error  in 
obtaining  failure  data,  (c)  using  finite  element  analysis  to  obtain  the  stress  field  throughout  a 
component,  and  (d)  using  numerical  reliability  algorithms  to  compute  component  probability  of 
failure.  The  discussion  here  focused  on  only  one  source  and  how  it  affects  component  proba¬ 
bility  of  failure.  It  is  not  the  intent  to  imply  that  this  type  of  loss  in  fidelity  is  paramount. 

The  rigorous  interrogation  of  the  C-ring  specimen  which  has  closed  fonn  stress  solutions  and 
effective  volume  and  area  solutions  available  has  helped  to  establish  an  analytical  protocol  for 
the  sectored  flexure  specimen.  Because  of  the  geometry  of  the  sectored  flexure  specimen,  closed 
form  solutions  are  not  available.  Thus,  the  procedure  established  for  benchmarking  the  C-ring 
specimen  was  applied  to  the  sectored  flexure  specimen.  Relative  to  the  sectored  flexure 
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specimen,  efforts  to  establish  an  asymptotic  value  for  ky-V  were  not  as  successful  as  with  ky -A . 
However,  establishing  a  value  for  ky'V may  not  be  as  important  as  the  k,t 'A  value.  Bend  bars  in 
general  tend  to  investigate  flaws  near  the  surface.  If  fractography  indicates  the  presence  of 
surface  flaws  along  the  outside  surface  of  the  sectored  flexure  specimens,  then  there  is  no  need  to 
fine  tune  the  computation  of  ky-V  for  this  specimen  geometry. 
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